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^^ . Abstract. — We derive the discrete version of the classical Helmholtz condition. Precisely, 

we state a theorem characterizing second order finite differences equations (see Definition [T| 
r^ I admitting a Lagrangian formulation. Moreover, in the afflrmative case, we provide the class of 

^\ . all possible Lagrangian formulations. 

C^ ' Keywords: Discrete Helmholtz condition; discrete calculus of variations. 



r~"^ ■ 1. Introduction 

> ■ 

^Q I A classical problem in Analysis is the well-known Helmholtz's inverse problem of the calculus 

»vq ' of variations (see [6l p. 71], [8] and [131 p.377]): find a necessary and sufficient condition 

under which a (system of) differential equation(s) can be written as an Euler-Lagrange 
^f^ ' equation and, in the affirmative case, find all the possible Lagrangian formulations. 

o' 

^^ I This problem has been studied by numerous authors and has been completely solved by 

A. Mayer [12] and A. Hirsch (91 110] . The formulation that we use is due to V. Volterra 

|14] . Precisely, let O be a second order differential operator. Then, the differential equation 

^ I 0{q) = can be written as a second order Euler-Lagrange equation if and only if all the 

H ■ Frechet derivatives of O are self-adjoint. This condition is usually called Helmholtz condition. 

- - -' We refer to [13] for a modern presentation and a complete proof of this theorem. 

A more difficult problem is to deduce from this characterization a complete classification of 
second order differential equations admitting a variational formulation. This has been only 
solved in dimension three by J. Douglas in his seminal paper [B] following a previous work of 
D.R. Davis [H [5]. We refer to [6l p. 74-75] and [HI p.377-379] for a historical survey. 

In recent years, an increasing activity has been devoted to a discrete version of the calculus of 
variations in the context of the geometric numerical integration. We refer to the book [Tj and 
the review paper [11] for more details. In this context, a second order discrete Euler-Lagrange 
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■ A_ (^{Q, -A+Q, T, h)) = 0, (1) 

for a given couple of Lagrangian (L_, L_|-) and where T is a bounded regular partition of M 
associated to the step size of discretization h. A_ (resp. A-|_) is the backward (resp. forward) 
finite differences operator associated to T. 

In this framework, we formulate the Helmholt'z inverse problem of the discrete calculus of 
variations as follows: find a necessary and sufficient condition under which a second order 
finite differences equation (see Definition [1]) can be written as a second order discrete Euler- 
Lagrange equation and, in the affirmative case, find all the possible Lagrangian formulations. 

This problem has been studied by numerous authors. We refer in particular to the work of 
Albu-Opris [1] and Cracium-Opris [3]. However, in each of these papers, the structure of 
the proof follows a different scheme than the continuous case and does not allow to make 
comparisons. Moreover, we enlarge these studies by taking account of non autonomous finite 
differences equations. 



2. Second order finite differences equations 

2.1. Partitions and finite differences operators. — In the whole paper, let us consider 
the following set: 

T:= {T = (tp)p=o,...,7V e M^+^ with TV > 4 and 3/i > 0, Vi = 0, ..., iV- 1, ti+i-ti = h}. (2) 

T is a set of bounded regular partitions T of M. Hence, for any partition T € T, an integer 
N = card(T) > 4 and a step size of discretization h > are associated. Consequently, for 
any T G T, we can also associate the following finite differences operators: 

A_ : M^+i -^ R^ (3) 



Q 

and 



h 



p=l,...,N 



A+ : M^+i -^ R^ (4) 



Q 



h 



p=0,...,N-l 

Let us note that A_ (resp. — A_|_) is the classical backward (resp. forward) Euler approxi- 
mation of the derivative operator d/dt. Moreover, A_ and — A+ commute and the discrete 
operator — A+ o A_ corresponds to the classical centered approximation of d'^/dt'^. 

Let us remark that all these previous discrete elements depend on T € T. For the reader's 
convenience, we omit this dependence in the notations. 
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In this paper, we are going to be interested in the discretization of differential equations 
defined on real intervals [a,b]. Hence, for any reals a < 6, we introduce the following set: 

Ta,b := {T G T with < to - a < h and < b - tN < h}. (5) 

For any reals a < b, Ta,b is then a set of regular partitions T of the interval [a, b]. 

2.2. Second order finite differences equations. — In the whole paper, let us note that 
we consider sufficiently smooth elements in order to make valid all the computations. The 
elements concerned are denoted by O, P, L, L„ and L+. 

In the continuous case, a second order differential equation on an interval [a, b] is defined by 
O"'' {q) = where O is a second order differential operator, i.e.: 

O: a<b ^ O"''' : '^^{[a,b],R) -^ 'r°([a,fe],]R) (6) 

q ^ O'^'^q) 

with: 

0"'*(g) : [a,b] ^ R _ (7) 

t ^ 0''^\qm = 0{q{t),q{t),m,t) 

where q (resp. q) is the first (resp. second) derivative of q and where: 

O: M^ — > R (8) 

{x,v,w,t) I — > 0{x,v,'w,t). 

Hence, a second order differential equation (independently of the interval [a,b]) is entirely 
determined by the application O. Let us give the following discrete analogous definition: 

Definition 1. — A second order finite differences equation, associated to a partition T G T, 
is defined by P'^{Q) = where P is a second order finite differences operator, i.e.: 

P : T € T ^ P^ : M^+i -^ R^-^ (9) 

Q ^ P^(Q) = (PJ(Q))^^,_^_, 

where 

\/p = 1, ..., iV - 1, Pj{Q) = P{Qp, (A_g)p, {-A+Q)p, (-A+ o A_g)p, tp, h) (10) 

and where 

P: R5 X M+ — > R (11) 

{x,v_,v+,w,t,^) I — > P{x,v-,v+,w,t,^). 
A second order finite differences equation (independently of the partition T & T) is then 
entirely determined by the application P. 

Let us consider a second order differential equation 0°''^{q) = on an interval [a,b]. Then, 
a usual algebraic way in order to provide a discretization of this equation is to consider a 
partition T G Ta^b and to define: 

P : M^ X R+ — > R (12) 

(a;,^;_,^;+,^«,i,,^) i — > 0(x, (1 — A)f_ + Aw^, w, t) 

with A G [0,1]. Hence, we obtain the numerical scheme P^{Q) = 0. The parameter A 
allows to choose for example the backward (A = 0), centered (A = 1/2) or forward (A = 1) 
approximation of the derivative d/dt. Such a discretization of 0"''''{q) = is said to be a 
direct discretization. 
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Example 1. — Let us consider the Newton's equation with friction q + q + q = defined 
on a real interval [a,b]. It is a second order differential equation associated to 0{x,v,w,t) = 
X + V + w. Hence, considering X = 1/2 and a partition T S Ta,b, we obtain by direct 
discretization the following numerical scheme: 



Qp+i — Qp~i Qp+i ~ '^Qp + ' 
~2h ^ 



Vp=l,...,iV-i, g^+ -^^+^^^-^^-^ + -^^+^ ^v,p^^P-i ^Q_ (^3) 



3. Formulation of the discrete version of the Helmholtz problem for second 

order finite diff"erences equations 

3.1. Reminder about the classical Helmholtz result for second order differential 
equations. — A continuous Lagrangian system derives from a variational principle. Pre- 
cisely, let us consider two reals a < b and the following Lagrangian functional: 

^«'* : 'tf^{[a,b],R) — > R (14) 

rb 

L{q,q,t) dt, 

where L is a Lagrangian, i.e. an application of the type: 

L: M^ — > m (15) 

{x,v,t) I — ;• L{x,v,t). 

Let '^Q{[a,b],M.) := {w G 'if'^{[a,b],M), w{a) = w{b) = 0} denote the set of variations. 
Then, q € ^^([a, 6],M) is said to be a critical point of ^"-^^ if for any variation w, we have 
D^°''^{q){w) = 0. A calculus of variations allows to characterize the critical points of ^"-'^ 
as the solutions on [a, b] of the following second order Euler- Lagrange equation: 

dL . , d f dL , ,\ ,^^r,h^ 

A dynamical system governed by such an Euler-Lagrange equation is called second order La- 
grangian system. We refer to [21 p. 55-57] for more details concerning continuous Lagrangian 
systems. 

The classical Helmholtz result is the following. Let O be a second order differential operator. 
The second order differential equation associated can be written as a second order Euler- 
Lagrange equation if and only if O satisfies the Helmholtz condition frequently given as the 
self-adjointness of all the Frechet derivatives of O"' for any reals a < b. We refer to |13j for 
more details. Nevertheless, the Helmholtz condition can be more explicitly formulated: O 
satisfies the Helmholtz condition if and only if 

Va<5, yqe'^\[a,b],R), ^ f^{q^q,q,t)\ =^{q,q,q,t). (H^ont) 

3.2. Second order discrete Euler-Lagrange equations. — Let us give the discrete 
analogous definitions and results of the previous section: 
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Definition 2. — A discrete Lagrangian functional, associated to a partition T G T, is de- 
fined by: 

N 7V-1 

Q ^ hY,L-{Qp,{A_Q)p,tp,h)+hY,L+{Qp,{-A+Q)p,tp,h), 

p=l p=0 

where {L^,L^) is a couple of Lagrangian, i.e. L± are applications of the type: 



L±: R^ X M+ 
{x,v,t,^) 



L±{x,v,t,^). 



(17) 



Let 



fN+l 



{W € M^^"^, Wo = Wn = 0} denote the set of discrete variations. Then, 
Q G ]R^+-^ is said to be a discrete critical point of S£^ if for any discrete variation W , we 
haveD^'^{Q){W) =0. 

A discrete calculus of variations allows to characterize the discrete critical points of ^"^: 

Theorem 3. — Let {L^,L^) be a couple of Lagrangian and T € T. Let ££^ be the discrete 
Lagrangian functional associated. Then, Q G M^"*"^ is a discrete critical point of ^^ if and 
only if Q is solution of the following second order discrete Euler- Lagrange equation: 

-^{Q, A_g, T, h) + -^{Q, -A+Q, T, h) 

+ A+ (^(Q, A_Q, T, h)\ - A_ (^(Q, -A+Q, T, h)\ = 0. (EL^) 

A discrete dynamical system governed by such a discrete Euler- Lagrange equation is called 
second order discrete Lagrangian system. 

Proof — Let Q G M^+i and W G M^+^ We have: 

TV 



p=i 



dL 



dL^ 



-Q^i*vWp + ^^{*p){A-W)p 



7V-1 



dL, 



dL, 



dx 



i**p)Wp + -^{**p){-A+W)p 



(18) 



where * := {Q, A_Q, T, h) and ** := {Q, —A+Q, T, h). Let us remind the following discrete 
integrations by part. For any {F, G) G M "'"^ x Rg ^ , we have: 

N N-l N-l N-l 

J2Fp{A,G)p=^{A+F)pGp and ^ Fp(A+G)p = ^ (A_F)pGp. (19) 

p=l p=l p=0 p=l 

Finally, combining (jlSp and (|19p . we obtain: 

AT-l 



p=i 



(9L_ , , dL+ . . ^ f dL^ . , , 
^(*.) + ^(**.)+A,(^(*)j^-A. 



which concludes the proof. 



^<"')J 


Wp, 


(20) 




D 



LOIC BOURDIN AND JACKY CRESSON 



Let US consider an Euler-Lagrange equation ( EL"' 1 defined on a real interval [a,b] and L 



the Lagrangian associated. Let us take for example L-{x,v,t,^) = L{x,v,t) and L+ = 0. 
Considering a partition T E Tab, we obtain that ^^ is a discrete version of ^"-^^ and 



( EL-^ ) is obtained by discrete variational principle on ^'^ . Such a method leads (JEL^J) to be 



called variational integrator: it is a numerical scheme for (jEL°'^ ) having the particularity of 
preserving its intrinsic Lagrangian structure at the discrete level. We refer to [Tl lllj for more 
details concerning the variational integrators. Let us note that one can also use a centered 
version by taking L-{x,v,t,^) = L+{x,v,t,^) = L{x,v, t)/2. 

Example 2. — Let us consider the Newton's equation without friction q + q = defined on a 
real interval [a, b]. It is a second order differential equation associated to 0{x ,v,w,t ) = x + w 



satisfying the continuous Helmholtz condition ( |Hcont| )- H corresponds to ( EL"'^ ) with the 
quadratic Lagrangian L{x,v,t) = (x^ — w^)/2. Considering a partition T G Ta.fe? taking 
L-(x,v,t,^) = L{x,v,t) and L+ = 0, we obtain the following discrete Euler-Lagrange equa- 
tion: 

yp=l,...,N-l, Q^ + 9p±1^1^1±9pz1 = o. (21) 

Let us note that ()2ip coincides with a direct discretization. Nevertheless, a direct discretization 
of an Euler-Lagrange equation do not lead necessary to a discrete Euler-Lagrange equation, 
see Example [^ In this case, we say that the Lagrangian structure is not preserved. 

3.3. Formulation of the discrete Helmholtz problem for second order finite dif- 
ferences equations. — Firstly, it is important to note that a second order discrete Euler- 
Lagrange equation is a second order finite differences equation (in the sense of Definition 
ED: 

Proposition 4- — Let {L-,L^) be a couple of Lagrangian. Then, the discrete Euler- 
Lagrange equation associated is a second order finite differences equation associated to: 

P : RS X M+ — > M. (22) 

{x,v-,v+,w,t,^) I — > ^ ^Aij^k{x,v-,t,^)v\ + Bij^kix,v+,t,Cy_ 



i+j+k>l 



where 



5L_ ^i+j+k-l Qi+J+k+lj^^ 

and where 

dL+ (-^y+j+k-i 9«+i+fc+iL+ 
B^,j,k{x,v,t,0 = 5(ij,fc)=(o,o,i)^^(^'^'*'0 + Tijj^i Q^^g^i+iQ^fc (^'^'^'g)' (24) 

where 6 is the Kronecker symbol. 

Proof. — We have just to take a partition T G T and to develop dL^/dv and dLj^/dv in 
power series. D 

We finally formulate the following discrete version of the Helmholtz problem: 

Discrete Helmholtz problem for second order finite differences equations: find a 
necessary and sufficient condition under which a second order finite differences equation can 
be written as a second order discrete Euler-Lagrange equation. Precisely, let P be a second 
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order finite differences operator. Our aim is to find a necessary and sufficient condition on 
P under which there exists a couple of Lagrangian (L„, L4.) such that for any T € T and any 
Q G M^+i, we have: 

P^iQ) = ^{Q, A_g, T, h) + -^{Q, -A+Q, T, h) 

+ A+ (^iQ,A.Q,T,h)] - A_ f^{Q,-A+Q,T,h)\ . (25) 



4. Solution of the discrete Helmholtz problem for second order finite differences 

equations 

Theorem 5. — Let P be a second order finite differences operator. Then, the second order 
finite differences equation associated can be written as a second order discrete Euler- Lagrange 
equation if and only if P satisfies the following discrete Helmholtz condition: 

(dP \ 8P 8P 

where * := {Q, A_Q, (-A+Q), (-A_ o A_Q),T, h) . 

Proof. — See Sections [5] and El D 



Let us note the similarity between the classical Helmholtz condition (Hcont ) and its discrete 



version ( Hjisc ) ■ It is also important to note that, similarly to the continuous case, the discrete 
Helmholtz condition ( HdJscD corresponds to the self-adjointness of all the Frechet derivatives 
of P^ for any T € T, see Section [S] 



Example 3. — Let us take the second order finite differences equation ()2ip obtained in Ex- 
ample\^ It is associated to the application P{x,V-,v+,w,t,S^) = x + w which satisfies the 



discrete Helmholtz condition ( Hdjsd ). It is expected because (j2ip is a discrete Euler- Lagrange 
equation by construction. 

Example 4- — Let us consider the differential equation q + sm{q)q = defined on a 
real interval [a, b] . It is a second order differential equation associated to the applica- 



tion 0{x,v,w,t) = X -\- sin{v)w satisfying ( [HcontD - It corresponds to (JEL"'''|) associated 



to the Lagrangian L{x,v,t) = {x'^/2) + cos(f). Let us consider T € ¥„;, and define 
P{x,V-,v+,w,t,^) = 0(2;, (1 — A)f_ + Xv+,w,t) with A G [0, 1]. Then, we obtain by direct 
discretization the following second order finite differences equation: 

Vp=l,...,iV-l,Q, + sm(^Q^+^ + (^-^^)^^^ + ^"-^)Q^-^)^^^^-\^/ + ^^-^=0. 

(26) 
Let P be the second order finite differences operator associated. Then, P does not satisfy 



( Hdisc ) O'^d consequently (|26|) can not be written as a second order discrete Euler- Lagrange 
equation. This is an example of direct discretization of an Euler- Lagrange equation not leading 
to a discrete Euler- Lagrange equation. The numerical scheme (|26|) does not preserve the 
Lagrangian structure of the differential equation at the discrete level. 
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5. Discrete Helmholtz condition and self-adjointness of Frechet derivatives of 

second order finite differences operator 

In this section, we prove that a second order differences operator P satisfies the discrete 
Helmholtz condition (Hdisc ) if and only if all the Frechet derivatives of P^ are self-adjoint for 



any T G T. 

5.1. Properties of the discrete derivative operators. — In this section, we first remind 
the classical discrete versions of the integration by part and the Leibniz formula. 

Lemma 6 (Discrete Leibniz formulas). — Let T € T and Q, W ^ M^"*"^. Then, we 
have: 

Vp = l,...,iV, (A_Qiy)p = (A_Q)pI^p + Qp_i(A_W)p (27) 

and 

Vp = 0,...,iV-l, (A+QW)p = (A+Q)pI^p + Qp+i(A+iy)p. (28) 

Finally, for any p=l,...,A^ — 1, we have: 

(-A+ o ^.QW)p = (-A+ o ^.Q)pWp + Qp{-A+ o A^W)p 

+ {-A+Q)p{-A+W)p + {A^Q)p{A^W)p. (29) 

For any T G T, let us denote R^,^^ := {W G M^+i, Wq = Wi = Wn-i = Wn = 0}. 

Lemma 7 (Discrete integrations by part). — Let T G T and (Q, W) G M^+^ x M^q^^ 

Then, we have: 

N-l N-2 N-1 N-2 

Y, QpiA-W)p = Y, {A+Q)pWp, and ^ Qp{A+W)p = ^ {A^Q)pWp. (30) 

p=l p=2 p=l p=2 

Finally, we have: 

N~l N~2 

J2 Qp{-^+ ° A_W)p = J2 (-^+ ° ^-Q)pWp. (31) 

p=l p=2 

5.2. Interpretation of the discrete Helmholtz condition as self-adjointness of 
Frechet derivatives of second order finite differences operator. — Let us define 
the following discrete version of the self-adjointness of a Frechet derivative for a differential 
operator: 

Definition 8. — Let P be a second order finite differences operator, T G T and Q G M^+-^ . 
We denote by DP'^{Q) the Frechet derivative of P'^ at the point Q. Finally, we denote by 
DP'^iQ)* the adjoint ofDP'^{Q) defined by: 

DP^iQY : E^+i -^ IR^-3 (32) 

Z ^ DP^{Qr{Z) = {DP^{Qr{Z))^^^_^_^ 

satisfying: 

N~l N~2 

V(l^, Z) G <o^i X IR^+\ h Y DP^{Q){W)Zp = hY, DP^{Qy{Z)Wp. (33) 

p=l p=2 

Finally, DP'^{Q) is said to be self-adjoint if for any p = 2, ...,N - 2, DP^{Q) = DP^iQ)*. 
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First, a simple calculation leads to the following result: 

Proposition 9. — Let P be a second order finite differences operator, T S T and Q € M ''~^. 
Then, for any p = 1, ..., N — 1 and any W G M^"'"^, we have: 



DP^{Q){W) 



dP 

dx 



dP 



dP 



dP 



i*p)Wp+-^i^p)iA^W)p+-g^{*p)i-A+W)p+—{*p)i-A+oA^W)p, 

(34) 
where * := (Q, A_Q, (-A+Q), (-A_ o A_Q), T, h) . 

Then, applying Lemmas [6] and [7] in Proposition [9l we can give explicitly the adjoint of a 
Frechet derivative: 



Proposition 10. — Let P he a second order finite differences operator, T G T and Q G 
M^+^ Then, for any p = 2, ...,N - 2 and anyWeR^^^, DP^{Q)*{W) is equal to: 



dP 

dx 



— (*p) 



dP 



dv^ 



A dP , , 



A — m^ - — 



+ 



p 
(*P-i) 



A dP , , 



(A_W) 



t^„ 



(35) 



+ 



dP 



dP 

dv- 



i*P+i) 



dP 



-A+W)p + — (*p)(-A+ o A_W)p, 



where • := (Q, A_Q, (-A+Q), (-A_ o A_Q), T, /i) . 

The main result of this Section is the following explicit characterization of second order finite 
differences operators P whose all Frechet derivatives are self-adjoint for any T G T: 

Theorem 11. — Let P be a second order finite differences operator. Then, DP^{Q) is self 
adjoint for any T G T and any Q G M^+-'^ if and only if P satisfies the discrete Helmholtz 



condition (Pdisc)- 



Proof. — Indeed, let T G T and Q G M^~^^. According to Proposition [TOl we have DP^{Q) 
is self adjoint if and only if the right term of (|34|) is equal to (|35p for any W G R^+^ and any 
p = 2, ...,N — 2. Consequently, DP^{Q) is self adjoint if and only for any p = 2, ...,N — 2, 
the three following equalities hold: 



A 9P, , 



p 
dP 



A dP 



dP 



A dP 

A^- — (* 



ow / ov 



A 9P^ , 



*P+i) 



dP 



(*P-i) 



dP 

dv+ 

dP 



dvA 



.*p)i 



0, 



di 



(*p), 



if and only if for any p = 2, ..., N — 1: 

dP 



A "w N 1 9P ^ ^ dP ^ 



The proof is completed. 



(36) 

(37) 
(38) 

(39) 
D 
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6. Proof of Theorem [5] 

6.1. Sufficient condition. — Let P be a second order finite differences operator associated 
to a second order discrete Euler-Lagrange equation. Let (L_,-L+) denote tlie associated 
couple of Lagrangian. According to Proposition U we have that P satisfies (j22|) . With a 
simple calculation, we can prove that for any T G T, any Q € M^"*"^ and any p = 1, ..., N — 1, 
the three following equalities hold: 






— ^((5p+i,(A_Q)p+i,ip+i,/i) - —^{Qp^i,{-A+Q)p_i,tp_i,h), (40) 



dP d'^L d'^L 



dv- dxdv p p / dxdv 



' 02 J 

A+{^-^iQ,A^Q,T,h)] (41) 



dP d^L d^L 

(*p) = ^Z^(<3p, (-A+Q)p,tp,/i) - ^-^((5p+i,(A_Q)p+i,tp+i,/i) 



dv4. dxdv dxdv 



A_(^(Q,-A+Q,T,/i)) , (42) 



where -k := (Q, A_Q, (— A+Q), (— A4. o A_Q), T, /i) . Finally, from these three equalities, we 
prove that P satisfies the discrete Helmholtz condition ( [HdiscD - 

6.2. Necessary condition. — Let P be a second order finite differences operator satisfying 



the discrete Helmholtz condition (H^isc)- The proof is based on the following proposition: 



Proposition 12. — Let P he a second order finite differences operator satisfying the discrete 



Helmholtz condition (Hdisc)- Let Li he the following augmented Lagrangian: 



Li : M^ X M+ — > R (43) 

{x,v^,vj^,w,t,^) I — > X I P{Xx, Xv-,Xv+, Xw,t,^) dX 

Jo 

and, for any T G T, let ^^ denote the following augmented Lagrangian functional 

^T . ^N+i _^ ^ (44) 

7V-1 

Q ^ /i^Li(Qp,(A_Q)p,(-A+Q)p,(-A+oA„Q)p,tp,/i). 
p=i 

Then, 

1. for anyT eT and any {Q, W) € M^+i x K^g^^ we have: 

N-2 

D^^{Q){W) = hY, Pp{Q)Wp, (45) 

p=2 
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2. there exists a couple of Lagrangian {L^,L^) such that for any T E T, any Q E R ^^, 
we have: 

Li(Q,A_Q,(-A+Q),(-A_ o A_Q),T,/i) 

= L_ (Q, A_Q, T, h) + L+ (Q, (- A+Q), T, /i) . (46) 

Proof — 1. Let T E T and (Q, VF) E M^+^ x ]RJ+^ We have: 

^i^(<5) = hy^Qp Ppi>^Q) dX. (47) 

Thus: 

N-l 1 AT-l 1 

L'^i^(Q)(W) = hS^Wp Pp{\Q) dX + hS^Qp DP^{XQ){XW) dX. (48) 



As P satisfies the discrete Helmholtz condition ( HdiscD and according to Theorem W\\ 
DP^{XQ) is self-adjoint. Using Definition [8l the fohowing equahty holds: 

7V-1 N-2 

h J2 DP^{XQ){XW)Qp = hY, XDP^{XQ){Q)Wp. (49) 

p=l p=2 

Then, from Equalities (j48|) and (|49p . we have 

N-2 1 N-2 1 

D^i^(Q)(Vr) = h^Wp PpiXQ) dX + hS^Wp XDP^{XQ){Q) dX (50) 

p=2 -^0 p=2 -^0 

Then, using d/dX[Pl'{XQ)) = DP'F{XQ){Q) and using an integration by part with respect 
to A on the second integral, we obtain: 

N-2 

DCj{Q){W) = hY, Pv{Q)Wp. (51) 

p=2 



2. Since P satisfies the discrete Helmholtz condition ( Hdjsc ) ; we have for any T E T and any 

Q E M^+i: 

1 fdP, , dP , A dP , , dP , 

where • := (Q, A_Q, (-A+Q), (-A_ o A^Q),T,h). As it is true for any Q E M^+^ we 
can differentiate the previous equality with respect to Qp-2 and Qp+i- It leads to the two 
following equalities holding for any T E T, any Q E R "'"^ and any p = 2, ..., A^ — 1: 

1 n d'^p . , 1 d^p. A 1 d'^p . . I d'^p . 

(*p-i) - To 7iZ2 (*p-i) = -T ^.. ^.. (*p-i) + T2 7^r:7^(*p-i) (53) 



/i \hdv-dw K^ dw"^ J hdv-dv+ K^ dwdv^ 

and 

I fi d'^p , ^ 1 a2p^ \ 1 52p ^ ^ 1 d^p , , 

(*p) + 72 tt::! (*p) = r7rr^r-(*p) + 7:2 7r77r-(*p)- (54) 



h V /i ai;+au; ' ^' ^ /i2 5u;2 ^'^^ / " h dv^dv. ' ^^ ^ h^ dwdv. '^ 
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Finally, we have for any T £ T, any Q G M ^-'^ and any p = 1, ...,N — 1: 

d^p , 1 f d^p ^ d^p , \ I d^p 

Hence, since Equality ()55p is true for any T € T, any Q S M "'"^ and any p = 1, ...,N — 1, we 
have for any {x,y, z,t,^) G M^ x M+: 

d^p f z-y \ I d'^p ( ^-y . c 



1 / d'^P ( z — y \ d'^P ( z — y \\ is 

Let us define: 

i: M^ X M+ — > M (57) 



2: — 11 

{x,y,z,t,^) I — > P[x,y,z,——,t,^ 



According to ()56p . we have: 



V(x,y,z,t,OeM^xM+, ^(x,y,z,t,O=0. (58) 

Consequently, the variables y and z are separable in i. Precisely, there exist two functions 

Q, /3 : R3 X M+ — >R such that: 

V(x, y, z, t, OeR^xRt, £{x, y, z, t, = a{x, y, t, + /3(x, z, t, 0- (59) 
Finally, we have for any T E T, any Q G M ~^^ and any p = 1, ...,N — 1: 

Li{*p) = Qp(fp{XQp,X{A.Q)p,\{-A+Q)p,X{-A+oA^Q)p,tp,h)dx\ (60) 

= Qp(f i{XQp,\{A^Q)p,X{-^+Q)p,tp,h)dx\ (61) 

= Qp(fa{XQp,X{A^Q)p,tp,h)+p{XQp,X{-A+Q)p,tp,h)dXj (62) 

= L_(Qp,(A_Q)p,tp,/i)+L+(Qp,(-A+Q)p,tp,/i), (63) 



where: 

L_ : M3 X M+ — ^ M (64) 

{x,v,t,^) I — > X / a(Ax, Au, t,^) dA 

JO 

L+ : M^ X M+ — > M (65) 

{x,v,t,^) I — > X / P{Xx, Xv,t,S,) dX. 
Jo 

The proof of Proposition [12] is now completed. D 



and 
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Now, from Proposition 1121 we are going to prove Theorem [5j Precisely, let us define, for any 
T G T, the following discrete Lagrangian functional: 

N N~l 

Q ^ /i^L_(Qp, (A_Q)p,tp,/i) +/i^L+(Qp,(-A+Q)p,tp,/i) 

p=l p=0 

where {L-,L^) is the couple of Lagrangian given in the point 2 of Proposition [T2j Then, we 
have for any T € T and any Q G M^+^: 



^' (Q) = ^l (Q) + /iL_ ( Qjv, ^" ^ '^^' ^ ) + ^^+ ( Qo' \ ^ *o, h\ , (67) 

where ^'^ is defined in Proposition [T2j Consequently, we have for any T G T, any Q G M ^^ 
and any W ^^^^^■. 

D^^{Q){W) = D^^{Q){W) = hJ2 Pp{Q)Wp. (68) 

p=2 

However, using the same method than in the proof of Theorem [3l we prove from Equation 
([66]) that for any T G T, any Q G M^+^ and any W G MJq^^: 

N-2 



p=2 






3x 9x \ dv J \ dv ^ 



Wp, 



(69) 
where * := (Q, A_Q,T,/i) and ** := {Q,-A+Q,T,h). Combining Equalities dMD and ([69]) . 
we conclude that for any T G T, any Q G M^'''^ and any p = 2, ..., N — 2: 

In order to finish the proof of Theorem O we have just to prove that Equality ()70p is still 
true for p = 1 and p = A^ — 1. We only prove it for p = N — 1. The case p = 1 can be proved 
in a similar way. 

In this way, let us take T G T and Q G M^^"*^. Let us denote a{T) = (tp+i)p=o,...,7V ^ T and 
(t(Q) = {Qp+i)p=o,...,N G M^''"^ where t^^i := tj\f + h and Qat+i := 0. From Equality ([70]) . 
we have: 

p;f-i(Q) = -pS(-w)) = ^(<'(*)«-2) + ^(-(**)«-2) 

where (t(*) := (o-(Q), A_o-(Q), o-(T), /i) and o-(**) := {a{Q),-A+a{Q),a{T),h) . Conse- 
quently, we have: 

P,L. («) = ^(-.) + ^("«-0 + A, (^,.)) ^_^ - A_ (^,„)) ^_^ . (72) 
The proof of Theorem [5] is completed. 
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7. Characterization of the null (couples of) Lagrangian 

In this section, we are interested in the second part of the Helmholtz problem both in the 
continuous and discrete cases. Precisely, once the Helmholtz condition satisfied, can we 
characterize all the possible (couple of) Lagrangian leading to the same second order (discrete) 
Euler-Lagrange equation? 

7.1. Reminder of the continuous case. — Let L^, L^ be two Lagrangian. They are 
said to be equivalent if they lead to the same second order Euler-Lagrange equation. In this 
case, we denote L^ ~ L^. The linearity of the Euler-Lagrange equation with respect to its 
associated Lagrangian implies that ~ defines an equivalence relation on the set of Lagrangian. 

Hence, the aim is to characterize the equivalence class of 0. If a Lagrangian L belongs to the 
equivalence class of 0, then it it leads to a null second order Euler-Lagrange equation in the 
sense that every curves q are solutions. In this case, L is said to be a null Lagrangian. We 
refer to |13| for a detailed proof of the following result: 

Theorem 13. — Let L he a Lagrangian. L is a null Lagrangian if and only if there exist 
two functions f : M? — > M and g : M — > M such that: 

Va <b,yqe <^2([a,6],M), L{q,q,t) = j^{f{Q,t)) +<?(*)• (73) 

Let us note that, in the major of literature, the previous theorem is presented with g = 0. 
Indeed, we have just to add an anti-derivative of g to /. However, in the next section, we can 
prove the discrete version of Theorem 1131 onlv with this presentation. 

7.2. The discrete case. — In the discrete case, we give the following discrete versions of 
the definitions and results of the previous section. 

Definition 14- — Let (LL,L^) and (L^,L^) be two couples of Lagrangian. We say that 
they a re equivalent if they lead to the same discrete second order Euler-Lagrange equation 
(EL^). In this case, we denote (LL,L^) ~ (L^,L^). The linearity of the discrete Euler- 



Lagrange equation with respect to its associated couple of Lagrangian implies that ~ defines 
an equivalence relation on the set of couple of Lagrangian. If a couple of Lagrangian {L-,L^) 
belongs to the equivalence class of 0, then it leads to a null second order discrete Euler- 
Lagrange equation in the sense that every discrete curves Q are solutions. In this case, 
{L-,L^) is said to be a null couple of Lagrangian. 

Our aim is then to characterize the set of the null couple of Lagrangian. It is done in the 
following result: 

Theorem 15. — Let (L_,L_|„) he a couple of Lagrangian. {L^,L^) is a null couple of La- 
grangian if and only if there exist two functions / : M^ x Wf — > M and g : W x Wf — > M 
such that for any T G T, any Q S M^"'"-'^ and any p = I, ..., A^.' 

L_ {Qp, (A_Q)p, tp, h) + L+ (Qp_.i, (-A+g)p_i, tp_i, /i) = A_ (/(Q, T, h))^ + g{tp, h). (74) 
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Proof. — Let us prove the sufficient condition. Let us assume that Equation ()74p is true and 
let ^ denote the discrete Lagrangian functional associated to {L-,L^) and to a partition 
T € T. Then, we have for any T G T and any Q G M^+i; 

N 7V-1 

^^(g) = hY,L-{Qp,{A^Q)p,tp,h)+hY,L+{Qp,{-A+Q)p,tp,h) 

p=l p=0 

N 

= h"^ ^L_{Qp,{A_Q)p,tp,h) +L+(Qp_i,(-A+Q)p_i,tp_i,/i) 
p=i 

= h^[A^{f{Q,T,h))^ + g{tp,h) 
p=i 

N 

= HQn, tN,h) - f{Qo, to,h) + h^ g{tp, h). 

p=i 

Consequently, since the set of discrete variations is Mq , every discrete curves Q G M "^"^ 
are discrete critical points of .if ^ and then the discrete Euler-Lagrange equation associated 
is null. Then, (L„, L+) is a null couple of Lagrangian. 

Now, let us prove the necessary condition. Let us assume that (L_,L+) is a null couple of 
Lagrangian. Then, for any T G T and any Q G M '''^, we have: 

BL BJ 

-^(Q,A_Q,T,h) + -^iQ,-A+Q,T,h) 

+ A+ (^{Q, A_g, T, h)) - A_ f ^(Q, - A+Q, T, h)) = 0. (75) 



Then, let us define: 



Bv J \ Bv 



M^ X IR+ — > m 

{xi,X2,t,^) I — > L_(xi, ^ ^ ,t,^ 



(76) 



and 

i+: ]R3 ^ M+ — > R (77) 

{xi,X2,t,^) I — > L+i^xi, ^ \ t,gj. 

Since Equality ([7?]) is true for any T G T and any Q G M^"*"^, we have for any {x, y, z, t, ^) G 



l^x 



* • 



— {x,y,t,0 + —{z,x,t + ^,0 + -^{x,z,t,0 + -^iy,x,t- CO = 0. (78) 

Then, by differentiating the previous equality with respect to y or to 2;, we obtain for any 

(x,y,t,0eM3xM+: 

{x,y,t,0+ ^ J" (j/,x,t-g,O = 0- (79) 



8x18x2 BxiBx2 
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Consequently, there exist two functions a, /? : M^ x 



p+. 



such that for any {x,y,t,^) G 



i.ix,y,t,0 + Uiy,x,t - to = ai^,t,C) + (3iy,t,C). (80) 

Moreover, let us denote Jf^ the discrete Lagrangian functional associated to {L-,L^) and 
to a partition T € T. Then, we have for any T € T and any Q € M^"*"^: 

N Af-1 

^^{Q) = hY^L^{Qp,{A^Q)p,tp,h)+hY,L+{Qp,{-A+Q)p,tp,h) 



p=i 

N 



p=0 



hY^ \L^{Qp,{A^Q)p,tp,h) +L+(Qp_i,(-A+Q)p_i,tp_i,/i) 

N 

h^ \i^{Qp,Qp_i,tp,h) +e+{Qp_i,Qp,tp-i,h) 
p=i 

N 



,^ \a{Qp,tp,h) + l3{Qp-i,tp,h) 



p=i 



With a discrete calculus of variations, we obtain for any T G T, any Q € M^"*^^ and any 



W € M^+i: 



N~l 



p=l 



D^''iQ){W) = Y, [Q^{Qp^tp^h) + ^{Qp,tp+,,h) 



dxi 



Wp. 



(81) 



Since (L-, L+) is a null couple of Lagrangian, we have D^'-^{Q){W) = for any T € T, any 
Q G M^+i and any W' G M^+\ Consequently, for any T G T and any Q G M^+\ we have: 



da '^(^ /r\ 

Tj — {Qp,tp,h,) + - — {Qp,tp+i,n) = 0. 

Consequently, we have for any (x,t,^) G M^ x M+: 

Hence, there exists a function 7 : M x R^~ — > M such that: 

V(x, t, C) G M^ X M+, a{x, t, + (3{x, t + C,0= lit, 0- 
Then, according to Equality ([5U|) . we have for any {x,y,t,0 G M'^ x M+: 

£- {x, y, t, + £+{y, x, t - ^, ^ = a{x, t, - a{y, t - C, + lit -^0- 
Thus, we have for any T G T, any Q G M "''^ and any p = 1, ..., A^: 

L_ [Qp, (A_Q)p, tp, h) + L+(Qp_i, (-A+(5)p_i, tp_i, h) 



(82) 

(83) 

(84) 
(85) 



a{Qp,tp,h) -a{Qp-i,tp-i,h) +j{tp_i,h) 



Hence, let us define: 



/ : ]R2 X 



(86) 



(87) 



ix,t,0 I — > ^a{x,t,0 
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and 

g: M X IR+ — > R (88) 

Then, we have for any T G T, any Q G M "^^ and any p = 1, ..., N: 

L„ {Qp, (A_Q)p, tp, h) + L+ (Qp_i, (-A+g)p_i, tp_i, /i) = A„ (/(Q, T, /i))^ + g{tp, h). (89) 
The proof is completed. D 
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